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Abstract 

We study the effects of a finite quark mass in the hedgehog configuration in the 
two phase chiral bag model. We discuss the chiral properties, such as the fractional 
baryon number and the chiral Casimir energy, by using the Debye expansion for 
the analytical calculation and the Strutinsky's smearing method for the numerical 
computation. It is shown that the fractional baryon number carried by massive 
quarks in the vacuum is canceled by that in the meson sector. A finite term of 
the chiral Casimir energy is obtained with subtraction of the logarithmic divergence 
term. 

1 Introduction 

Strangelets of finite volume quark matter with strangeness has been one of the most 
interesting subjects of exotic particles in hadron and quark physics ^12113111. In recent 
studies of strangelets, a model has been presented for a discussion of chiral symmetry 
breaking inside a quark droplet of finite volume quark matter [HI El 13 E]- There, in 
addition to quark confinement for a finite volume system, it was allowed that the quarks 
acquire a finite mass by dynamical chiral symmetry breaking. In this picture, it was 
assumed that the quark interaction was provided by the Nambu-Jona-Lasinio (NJL) type 
interaction inside a quark droplet. In earlier works, the quark wave function was given by 
the MIT bag in Refs. [3 El EI- Later, in order to overcome the chiral symmetry breaking 
at the bag surface, the chiral bag model was introduced in Ref. . The latter model was 
called as the NJL chiral bag model. Such an idea was first presented by T. Kunihiro in 

1983 pun]. 

In the discussion of strangelets by the NJL chiral bag model, the hedgehog ansatz 
was used for the pion and ud quarks to include non-linear interaction of pions and quarks 
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|H]. In the history of the study of the chiral bag model since [HI [12], effects of vacuum 
polarization, such as the anomalous baryon number and the chiral Casimir effects, have 
been discussed extensively for massless quarks [13 [3 [13 [IHl [HI [HI • It was shown that 
the correct behavior of the baryon number was provided by the contribution from the sea 
quarks [13 [THl [SOI- In the development of the chiral bag model, several techniques for 
calculation of the sea quark contribution was also developed. In Ref. ^Zj, the fractional 
baryon number was derived by using the Multiple Reflection Expansion method. After 
that, the Debye expansion was used P3[2nj- The Casimir energy is also affected by the 
pion cloud. It was shown that the chiral Casimir energy had a term of logarithmic diver- 
gence [221[2ni- However, these properties for massive quarks have not been investigated 
so far. In the framework of the NJL chiral bag model, the dynamical quark mass is an 
order parameter for the chiral symmetry breaking in the chiral bag. In order to search 
the real stable state in the energy variation with respect to the dynamical quark mass, we 
need to clarify the effect of vacuum polarization for the finite quark mass in the hedgehog 
configuration in the chiral bag. This is the motivation in the present paper. 

We organize this paper as follows. In Section 2, we propose a model lagrangian of 
the chiral bag with massive quarks. In this paper, we describe the quarks with a finite 
Dirac mass, which is generated through the NJL type interaction in the chiral bag. In 
Section 3, we discuss a fractional baryon number for massive quarks by using both of 
numerical computation and analytical formulation. In Section 4, we discuss the chiral 
Casimir energy. Deriving the terms of logarithmic divergence the finite term of the chiral 
Casimir energy is obtained numerically. We summarize our discussion in Section 5. 

2 The chiral bag with massive quarks 

In the NJL chiral bag model [HI [HIE], the quark mass in a finite volume system is induced 
by the NJL point-like interaction |2S|- In the NJL chiral bag model, we consider the ud 
quark sector with a finite quark mass induced by the mean field approximation in the 
scalar channel. The lagrangian in the ud quark sector is written as 



which includes a boundary condition explicitly. Here, ip = {u, dY is the ud quark fields, 
and the ud quark mass matrix is given by m = dmg{mu,m(i) with flavor symmetry 
m = rriu = irtd. The ud quark mass is the sum of the current mass and the constituent 
mass which is generated in the non-zero expectation value of the quark scalar condensate 
qq. In general, the constituent quark mass can be position dependent. Here, in order to 
simplify the essential discussion of the finite quark mass, we treat the quark mass as a 
constant. For a bag configuration, we assume a static spherical bag with radius R. The 
step function is multiplied in the first term in order to confine quarks inside the bag. 
Here, r is a distance from the center of the bag. The second term with the delta function 
realizes the chiral invariant boundary condition at the bag surface, where we define 
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with the TT meson field vr and the Pauh matrix r. In this paper, we do not exphcitly show 
the lagrangian in the pion sector, since our current interest is to study the chiral vacuum 
polarization of the quarks in the bag. 

In order to consider the non-linear effect of the pion, we assume the hedgehog ansatz 
in the tt meson sector, where the pion field conserves the grand spin K = J + 1 with total 
angular momentum J and isospin I. The hedgehog pion field for r > R is written as 



TT 



F{r)n, (3) 



with a chiral angle F{r) and a unit radial vector n in the real space 

According to the n meson sector, we introduce the hedgehog basis set in the ud quark 
sector [iHl EZI- We construct the quark wave function ■i/'^"^ for natural (k = +) and 
unnatural (k = — ) assignment for grand spin K, respectively; 



^(+) ^ f aojK{pr)\0) + aijK{pr)\l) 

I a2jK+iipr)\2) + a3jK-iipr)\3) 



(4) 



and 



,/,(-) _ ( b2jK+i{pr)\2) + b3jK-i{pr)\3) \ , . 

^ 1^ bojKiprm + hjK{pr)\l) )■ 

The coefficients aj and 6, (z = 0, ■ ■ ■ , 3) are determined by satisfying the equation of 
motion of the ud quark for k = ±, respectively, jxip^) is the spherical Bessel function 
and p is the ud quark momentum. Here, the two-component spinors |0), ■ ■ ■ , |3) are given 
by 

|o) = YKMie,^)xl 

|1) = {KM - ixli^\KM)YKM-,ie,<P)xl 

Ai=-l,0,l 

|2) = {K + 1M - filfi\KM)YK+iM-f.{e,<P)xl 

Ai=-l,0,l 

|3) = J2 {K -IM - f,lf,\KM)YK^,M-,{eA)xl (6) 

Ai=-l,0,l 

where YiMi^, <P) is the spherical harmonics with spherical coordinate (6', 0). are eigen- 
states oiG = S + T=a/2 + r/2, 

xi = ^{\m)-\i)\u)). 

X\ = \])\u), 

Xl = ^{\m) + \i)\n)), 

X\ = \i)\d). (7) 
The sign of naturalness k corresponds to the parity P = {—)^^'^. 
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In order to satisfy the equation of motion, the coefficients Oj for naturalness assignment 
K = + are subjected to 



( E — m 


iPp 



—iPp —iQp \ 

E — m —iQp iPp 

iQp E + m 

-iPp E + m J 



f o,o\ 

V as y 



0, 



(8) 



where we define 



P 



Q 



I K+1 
2K + 1' 

I K 



E — ^Jp"^ + w?. 
Then, we obtain two independent solutions 



(9) 



a' = 



N' 



( ^iP \ 

E—m 

^iQ 

E—m ^ 
1 



V 







N" 



E—m 



E-m 





iQ \ 
^iP 



(10) 



/ 



Here, the normalization constants N' and N" are determined by the condition / d^xip^'^^'^tlj^'^^ 
1, 



E 



E-m 

E 
E-m 



jKipRf + 3K+,{pRf - ^ (2(X + 1) - I) 3K{pR)3K+i{pR) 
jK-iipRf + Jk{pR? -^(2K+ I) jK-i{pR)jK{pR) 
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The final solution is expressed as a linear combination of a' and a", 



a = c'a' + c"a" = c'N' 



( -^iP \ 

E—m 

-^tQ 



E—m 
1 





-^ip 



J 



E—m 





V 



:i2) 



with constants d and c" . We obtain a solution for the unnaturalness assignment k — 
in the same way. 
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The eigenenergy of the hedgehog ud quark is given by the boundary condition at the 
bag surface r = R, 

in ■ 7V;(") = -e'^^-"^(^)^5V^("). (13) 

By substituting the solution of Eq. (fT^ . we obtain an equation for the eigenvalue, which 
is given by 

cosF{R)!^jK{pRf- jK+i{pR)]K-i{pR) 

-/« — - — jKipR) {jK+l{pR) -jK-i{pR)} 
E — nm sin F(R) , ^ r , ^ / ^ , ^ 

+ ^^l/ ji^ {jK+iipR) + jK-iipR)} = 0. 14 

p 2K + I 

The equation for = is obtained by setting = 0. Eq. (fT^ determines the 

energy spectrum of a quark in the bag. We mention that there is a symmetry of the 
energy level in the positive and negative energy 

Ekp{F) = -Ek-p{-F). (15) 

This relation has its origin in the invariance of the lagrangian (P) under the transformation 
U ^ U* or F —>■ —F. Another symmetry under F — > F + vr and k — > — k is conserved in 
the massless case. However, this symmetry is not conserved for the finite mass. This is 
due to an asymmetry of the quark energy levels, which is shown in Fig. 1 in [Hj. 



3 Baryon number conservation 

In the limit of the zero bag radius, the meson field in the Skyrme lagrangian gives a 
mapping from R^ ~ to SU (2) ~ by imposing a boundary condition [/ = 1 at 
r ^ oo. The winding number in the mapping is identified with the baryon number 
[2H1 12nj • The fractional baryon number due to the pion cloud is given as |2Sl IIH] 

1 



-Btt = 

TT 



sin 2F 
2 



(16) 



where F is the chiral angle at the bag surface. In the pure Skyrmion, the chiral angle is 
given as ravr with integer n. 

On the other hand, in the chiral bag with finite bag radius, the chiral angle is not 
generally equal to rni. Therefore, the meson carries only a fractional baryon number. 
There, it has been known that the total baryon number is composed of the fractional 
baryon numbers of the pion, of the vacuum and of the valence quarks in a bag fHl I2Z] • 
The baryon number carried by the vacuum quarks is defined as 

5,(m,F) = i(0|[^t,^]|o) 

= -^lim^sgn(F„)e-*l^"l (17) 
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Here, \0) is a vacuum state of an empty bag filled with negative energy quarks, En = 
En{m, F) an n-th state energy of a quark with mass m, and at the chiral angle F at the 
bag surface. The quantum number n labels the grand spin and parity. The sum is taken 
over for all the quark states n with positive and negative energies in the bag. In order to 
obtain the convergence in the sum, an exponential type regularization is multiplied. 

In the studies of the chiral bag model with massless quarks, it was shown that the 
baryon number of the vacuum quarks was canceled by that of the meson cloud fSJ ^3 
113 HZl UHl 120] • Accordingly, the total baryon number is a conserved quantity. This 
result should also be the case for the chiral bag with massive quarks. It is expected that 
the baryon number conservation is not affected by the finite quark mass, because U{1)b 
symmetry is conserved in our lagrangian (^. In the following, we show explicitly the 
cancellation of the fractional baryon number between the vacuum massive quarks and the 
pion cloud. We use both numerical computation and analytical procedure. 

First, we show the numerical computation. In the studies of chiral bag model with 
massless quarks, the regularization, such as the Gaussian type j22] , heat kernel type j2E] , 
have been used with success. In this paper, we use the Strutinsky's smearing method 
|24[ 123] . This method has an advantage that the states necessary for computation are 
limited up to the grand spin K^ax ~ 40 in order to obtain a good convergence, while we 
need K^ax ~ 100 at least for the other regulators. We rewrite the baryon number (fT7|) 
by introducing a delta function for a density of states to pick up the discrete levels 



1 r°° 

Bg{m, F) = -- dxY, sgn(K)5(x - E„i?) 

2j J— oo „ 



Then, we replace the delta function by a gaussian function, 

1 

Bg{m,F) = -- dxp,{m,F), (19) 
where we define a function for a density of states 

p.{m,F) = ^sign(K) y= , (20) 



where 7 ~ 2 — 3 is a smearing parameter. We obtain Bg = at F = 0, since the energy 
spectrum is symmetric for positive and negative energies. In oder to obtain a rapid 
convergence in the sum, it is convenient to define the difference of Px{iTi, F) between finite 
F and F = 0, 

Px{m,F) = p^{m,F) - p^{m,0). (21) 
Therefore, the baryon number defined by Eq. (fTTjl is given by an alternative formulation 

/oo 
dxp,{m,F). (22) 
-00 
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For a numerical computation, we restrict the range of the x-integral in a finite interval 

^ [ ■^maxi ■^max\- 



B,im,F) 



dxpxim, F). 



(23) 



We use Xmax — 20 to obtain sufficient convergence. Carrying out numerical computation, 
we obtain the result which is given by 



B,im,F) = - 

71 



F- - sin2F 
2 



(24) 



for F G [Fq, Fq + n]. The value of Fq depends on the quark mass; Fq/ti = —0.5, —0.695, 
—0.816 and —0.878 for mR = 0, 1, 2 and 3, respectively. The equation is independent 
of quark mass m, and coincides with the opposite sign of Eq. |T6|l in the meson sector. 
Therefore, the fractional baryon number by meson and vacuum quark cancel each other 
for any quark mass. 

The total baryon number is equal to A for a quark droplet with 3A valence quarks for 
any chiral angle. Let us investigate the case of a baryon. The baryon number carried by 
valence and vacuum quarks for chiral angle Fq < F < is given as a sum of the valence 
quark and the fractional baryon number 



B{m, F)yal+sea = 1 H 

TT 



F - - sin 2F 
2 



(25) 



For —71 < F < Fq, the valence quark is absent since the 0^ state is absorbed in the 
vacuum. Then, the baryon number is obtained by substituting F + tt instead of F from 
a periodicity. 



B{m, F)^ai+s 



1 

vr 

1 + 



(F + tt) - -sin2(F + 7r) 



vr 



F - - sin 2F 
2 



(26) 



Therefore, the total baryon number is one for any chiral angle. In the same way, the 
conservation of baryon number holds for any quark droplets. 

Second, we present an analytical calculation of the baryon number ()17|1 by using the 
Debye expansion. This method was originally developed in the chiral bag with massless 
quark ED] • We apply this formalism to the massive case. We rewrite Eq. (fTTjl in terms 
of the quark propagator as in [13 120] • Then, we obtain 



B,{m,F) 



- lim 

2 r^0+ 



— e*™ y (2ir + l)^ln 



K=0 



dx 



Sxijn, F; ix) 
S*j^{m, F; ix) 



(27) 



Here, Sxijn, F;p) is given as a product of the left hand side of the eigenvalue equation 
Eq. (fT^ for natural {n = +1) and unnatural (k = —1) states. For the sake of using the 



7 



Debye expansion, we substitute an imaginary number ix for pR in the spherical Bessel 
function. We exphcitly show 

So{F;zx) = cosF [Jl{^xf - U^xf) - 2 Uix)ji{ix) (28) 



P 



for K = 0, and 



Sxijn^ F; ix) 



— cos F (jK{ixf - jK+i{ix)jK-i{ix)) + J^r^^J xiix) {jK+i{ix) +jK-iiix)) 



p 



+ 



m 



P 



- 2 

{jKiixf + jK+i{'i'X)jK-i{ix)^ -Jk {jK+i{ix) -jK-i{ix)) 



for K ^ 0. For a short notation, we write rh = mR, p = x and E = \/fa^ + x^. The 
x-integral picks up an eigenvalue as a residue instead of solving directly the boundary 
condition. We note that Eq. (fTTj) is an exact formulation, but it is not practical to 
consider the x-integral and the sum over K without approximation. In order to obtain 
the final result, it is sufficient to use an asymptotic behavior of the quark energy spectrum. 

In the following manipulation, we discuss the K = ^ and K \ components separately 
for convenience. First, we consider the = component, which is defined by 



1 



Bf='^\m,F) = - lim 



2 r^04 



" — e^™-^ln 
-oo 2i7i ax 



So{m, F; ix) 
S'o(m, F; ix) 



It is convenient to write the function S'o(m, F; ix) in a polar coordinate 

So{m,F-ix) = i?o(a;)e**«(^). 
By using the asymptotic form 

ji{ix) ~ ^(e'^ + e" 

we obtain 



arctan 



^2x 



e ^/ x^ + m? sin F + m 



e^^ + 2 



-2x 



xcosF 



Then, it is straightforward to perform the integral, giving 

i?f=°)(m,F)^^^, 

TT 



(30) 



(31) 



(32) 



(33) 



(34) 



for any quark mass m. 

Second, for K > 1, we use the Debye expansion for the modified Bessel function in 
estimation of Eq. ()29|1 . In our notation, the modified Bessel function Iu{x) is defined by 



JK[IX) 



-fUx) 



(35) 



(29) 
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with V = K + 1/2. The Debye expansion is a uniform asymptotic expansion with no 
constraint between u and x. The Debye expansion gives 



7TV2t 



(36) 



in the lowest order, where we define t = and f^{x 

Now, let us define 



5f = - lim r — e*"" y 



E'>i 



2 r^0+ 

Here, we write Sk{iti, F; ix) in a polar coordinate 



t — u sinh ^(z//x). 

Sxifn, F; ix) 
S'^{m,F]ix) 



(37) 



(38) 



Then, we apply the Debye expansion (jHBj) . After a little tedious calculation, we pick up 
only the terms of the leading order of m, and obtain 



arctan 



2K+1 V2t3 



^sin2F 



4z/m\ 



(39) 



We consider that the quark mass m is smaller as compared with the momentum p. Indeed, 
this is a good approximation, since we are interested in the asymptotic behavior of the 
quark energy spectrum. The x-integral for ^k{x) is integrated out and the integral is 
determined only by ^k{oo) — ^k{—oo). Concerning the term which has no quark mass, 
we use an identity 



^ cos[Tiyz) _ _i^^^^(^^/4) _ 2cos(rz/2), 

K>1 ^ 



and 



lim 



°° dx 
-oo 27ri 



K>1 



dx 



—I 



ux 
2t3 



1 

'2^' 



(40) 



(41) 



Concerning the term proportional to m, we use 



lim 



oo 27ri dx t 



(42) 



K>1 



which indicates the mass term does not contribute to the integral. Consequently, we 
obtain the result 



Bf^^\m,F) 



sin2F 
2ti 



(43) 



Finally, by Eqs. ()34p and ()43|) . the fractional baryon number obtained in the analytical 
procedure coincides with the numerical result Eq. ()24|1 . 
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4 Chiral Casimir energy 

The chiral Casimir energy arises as a result of the modification of quark energy levels for 
finite pion cloud. Following the regularization scheme in the baryon number, the chiral 
Casimir energy is defined by 

Ecim^F) = ^{0\[^P\H^]\0) 

= limEsign(i^n)i^^ne-*l^"l, (44) 

as in Refs. CI 113 dl HH 123 123 • 

We show the numerical procedure by using the Strutinsky's smearing method. By 
using the Gaussian function (j^ . we rewrite Eq. (jl^ as 

/oo 
dxxpx{m,F). (45) 
-oo 

By choosing the reference point at F = 0, we obtain an alternative formulation 

AEc{m,F) = Ec{m,F)-Ec{m,0) 

/oo 
dxxpx{rrt,F). (46) 
-oo 

Restricting the range of integral in a finite range of [—Xmax^Xmax] with sufficient conver- 
gence, we obtain, 

dxxp^{m,F). (47) 

In order to achieve a rapid convergence of the integral, it is much more practical to use 
the partial integration^. We mention that we need Xmax ~ 40 to obtain the convergence 
in the calculation. More states are necessary for the chiral Casimir energy as compared 
for the baryon number {x^ax ~ 20). 

For massless quarks, it has been known that there is a logarithmic divergence term 
in the Casimir energy (03}. It was shown that the divergence term is proportional to 
sin^ F [ini El Uni J-9, 20^ EH 1201 12Z|- Here, we also derive the logarithmic divergence 
term for massive quark by using the Debye expansion developed in |19| 120] • Following the 
procedure in ^31201; we rewrite the Casimir energy as 

1 r°° dr °° d 

AEcim, F) = — hm / ^xe^^ ^ {2K + 1);^ In 

Just as in the case of the baryon number, we consider the calculation for K = and 
K > 1, separately for convenience. 

^We use a relationship dxxpx{rn, F) = — dx^p'^{m, F) with Pxim, F) ^ for |a;| ±oo. 



Sxim, F; ix) 



Sxim, 0; ix) 



(4J 
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First, we investigate the K = sector, 



AE. 



{K=0) 
C ' 



1 



m,F) = — lim / — a;e*™— In 



2R r-O 



dx 



d 



2tx 



dx 



So{m, F; ix) 



Soim, 0; ix) 



(49) 



We apply the asymptotic form of the spherical Bessel function ()H2j) . By expanding 
\So{m, F; ix)/So{m, 0; ix)\ in the lowest order of m, we obtain 



1 



m, F) ~ — lim / — xe*™ — 



2R r^o+ 



°o dx 



d 



271 



dx 



2777, 

cos^F + sinFtanh^(2x) +sinF — tanh^(2x) .(50) 

X J 



The first two terms in the integrand, which does not contain quark mass, are given in 
[T^ 120] . The integral of the third term proportional to rh is easily shown to be equal to 
zero. Consequently, we obtain the result 



AFjf-°^(m,F) 



AttR 



(51) 



The divergence term does not appear in the K = sector. 
Second, we investigate the sum over K > 1, 



AE, 



c ' 



m, F) 



— lim 

2R r-0+ 



°° dr °° d 

— xe*™ y (2ir + l)-^ln 

-oo 271 dx 



SKim, F; ix) 



Sxim, 0; ix) 



(52) 



After a little troublesome calculation by using the Debye expansion and expanding in the 
lowest order of m, we obtain 



In 



Sxifn, F; ix) 



SxiiTi, 0; ix) 



~ — sin^ F 

,2 



+ 



1 



l + Ll){l + L^_) (2K + 1)2(1 + L^)(1 + L? 



m 



2x2 



sin^ F- 



;i-L+L. 



1 + L,)2(1 + L_ 



— sin^ F— 
2 



+ — sin^F- 1 



t2 / t2 



t 



(53) 



Here, we define L± = for a short notation. The first term, which does not 

contain quark mass, coincides with the case of massless quarks fOlEHI- The second term 
proportional to fh is further analyzed as followings. By taking a derivative with respect 
to X, we obtain 



^ 1 
X— In 

dx 



Sxim, F; ix) 



Sxim, 0; ix) 



1 2 ^ / x^ 6 31z/2 



+ - 



— 2 

m 



sin^F 



2 



3 



1^ 



46z/^ 21z/6 
+ 



t7 

23z^ 



t9 



41z/^ 21z/f\, 



t7 
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Then, introducing new variables h = t/v and z = x/v, we have 

(i^>i) sin^F /-oo ^ 

A£/>. - '(■m,r) ~ -— iim / dz—^ > coslrz/^l 

^ ^ ^ vri? T^o+Jo 

sin^F , f Q 31 46 21\ ^ cos 



• 2 r. 1- /"^ , / 1 6 4 \ ^ , , 

+— — sm F hm / rf^; -— + 7-7-7^ > cos(rz/2;) 

m2 r , f ^ 23 41 21\ ^ cos{tvz)^. 




The summation over K > 1 and the integral is achieved by the formula Eq. (j40|l . Finally, 
we arrive at the result 

,^(K>i). r 3 / 2 , 121n2-9\ 

Therefore, we obtain a logarithmic divergence term with coefficient of sin^ F for a finite 
quark mass. We note that our result reproduces the previous result for the massless 
quarks in m — >^ [HI EOj • 

Now, in order to remove the logarithmic divergence, we subtract the second derivative 
of the chiral Casimir energy with respect to the chiral angle at the reference point F = 0, 
and we obtain a finite contribution 



Ei:''{m, F) = AEc{m, F) - - sin' F 



2 dF^ ^=0 



(57) 



In Fig. ^ we show the chiral Casimir energy F^*" as a function of the chiral angle F. 
The lines are distinguished by the quark masses mR = (dashed line), mR = 1 (solid 
line), respectively. It is a remarkable point that the Casimir energy for the massive quark 
takes a nonzero value at F = — vr, while that of the massless quark becomes zero at 
F = —71. This is because that the massive quark has an asymmetric energy spectrum 
between F = and — vr, while the massless quark has a symmetric spectrum, as shown in 
Fig. 1 in 0. 

Here, we recall the Cheshire Cat picture in the chiral bag model with massless quarks. 
There, the continuous transformation from the chiral bag to the Skyrmion was induced 
in the limit of small bag radius ^Sl ^1 ^j- This picture seems not to be applied to 
the case of the massive quarks, since the chiral Casimir energy takes a finite value at 
the zero bag radius, or F = — vr. However, this observation does not make us abandon 
our discussion. Indeed, in our formulation, the quark mass is a dynamical variable which 
should be determined by the energy variation. Our previous result suggested that the 
chiral symmetry was restored in the limit of small bag radius |HI. Therefore, the NJL 
chiral bag model holds the Cheshire Cat picture. 



12 



-2 



-1 

F/7I 







Figure 1: The Casimir energy as a function of the chiral angle at the bag surface. The dashed 
and sohd hues indicate the quark masses mR = and mR = 1, respectively. 

5 Summary 

We discussed the chiral properties in the chiral bag model which contained massive quarks. 
Considering that the dynamical quark mass was generated by chiral symmetry breaking, 
the fractional baryon number and the chiral Casimir energy were investigated. We dis- 
cussed the effects of the finite quark mass in the hedgehog configuration with an emphasis 
on the technical procedure. 

We showed the numerical calculation by the Strutinsky's smearing method and the 
analytical technique by the Debye expansion. It was shown numerically and analytically 
that the fractional baryon number carried by vacuum massive quarks inside the bag was 
canceled that of the vr meson outside the bag. Therefore, the total baryon number is 
exactly conserved. By using the Debye expansion, it was shown that the chiral Casimir 
energy had a logarithmic divergence term. The chiral Casimir energy was obtained nu- 
merically by the Strutinsky's smearing method by removing the divergence term. It is a 
point different from massless quark that the chiral Casimir energy for massive quark had 
a nonzero value at the chiral angle F = — tt. 

As further development, we plan to discuss a fully chiral symmetric equation of motion 
for a finite quark mass. In our present analysis of the NJL chiral bag model, we performed 
a mean field approximation only in the scalar channel in the four point interaction jS]. 
However, in the hedgehog configuration, it may be allowed to have a finite expectation 
value of qq, not only in the sigma channel, but also in the pion channel. In the latter 
case, the expectation value of qq with pion quantum number can be given as a finite value 
in a basis set of the hedgehog quark wave function. This study is now in progress. The 
quantization of the hedgehog configuration for massive quark is an important subject in 
order to obtain the state with definite spin and isospin. It is also an interesting subject to 
include the strangeness sector in our framework. These subjects are left as future works. 
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